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Abstract: 

Using new experimental data on the leptonic mass spectrum of B — > niv, we simultaneously 
determine \V u b\ and constrain and aj, the first two Gegenbauer moments of the pion's 
leading-twist distribution amplitude. We find \V u b\ = (3.2 ± 0.1 ± 0.1 ± 0.3) x 10~ 3 , where 
the first error is experimental, the second comes from the shape of the form factor in q 2 
and the third is a 8% uncertainty from the normalisation of the form factor. We also find 
af(l GeV) = 0.19 ± 0.19 and aj(l GeV) > -0.07. 
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\V u b\ is one of the least well-known elements of the CKM quark-mixing matrix. A more 
precise determination of this parameter will not only greatly improve the constraints on 
the unitarity triangle, but also provide a stringent test of the CKM mechanism of flavour 
structure and CP violation. In this letter, we determine \ V u b\ from the exclusive semileptonic 
decay B — > tt£u, based on the invariant lepton-mass spectrum recently reported by BaBar 

and the light-cone sum rule calculations of the relevant form factor in Ref. [2|. 

The hadronic matrix element relevant for B — > niu is given by 

(n(p^\m,b\B{p B )) = (vB+P.-q^^J U^ + ^^q.Hq 2 ), (1) 

where the form factors / +i o depend on q 2 = (ps — Pn) 2 , the invariant mass of the lepton-pair, 
with GeV 2 < q 2 < (tub — m n ) 2 = 26.4 GeV 2 . /+ is the dominant form factor, i.e. the only 
one needed for calculating the spectrum in q 2 , 

rlV C 2 IT/ I 2 

^-^^W^Mf)? (2) 
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for massless leptons; \(q 2 ) = {rn 2 B + m 2 — q 2 ) 2 — 4m B m 2 is the usual phase-space factor. The 
determination of \V U \,\ from B — > -niv requires theoretical input on / + , which has been the 
subject of many a calculation using various methods, in particular quark models PJ, QCD 
sum rules on the light-cone (LCSRs) IH E] and lattice simulations [6J- The challenge for 
theory is twofold: the region of applicability of theoretical calculations is, in most cases, 
restricted to part of the full physical phase-space; a calculation of /+ for, say, small values 
of q 2 is, however, not sufficient, as experimental data on the decay spectrum are still very 
scarce, so that any meaningful extraction of \V u b\ necessitates the extrapolation of the form 
factor to all q 2 . LCSR calculations, for instance, are valid for large pion momentum, which 
translates into small to moderate q 2 < 14 GeV 2 , whereas lattice calculations are restricted to 
small pion momentum, corresponding to large q 2 £ 15 GeV 2 . Extrapolations rely either on 
a model for the g 2 -dependence of /+, like vector meson dominance or the parametrisation 
advocated by Becirevic and Kaidalov [7j, or dispersive bounds on the form factor, which have 
been studied for instance in Ref. [Hj- The BaBar collaboration has measured the spectral 
decay distribution in 5 bins in q 2 |Tj, which is a significant improvement over previous results 
reported for 3 bins 9], and allows one, for the first time, to assess the validity of various 
parametrisations of the g 2 -dependence of /+ like 

• vector meson dominance (VMD); 

• the parametrisation of Becirevic and Kaidalov (BK) [7j; 

• the extended BK parametrisation used by Ball and Zwicky (BZ) [2]. 

All these parametrisations can be motivated from the exact representation of / + in terms of 
a dispersion relation, 

Res g2=m | J+{q 2 ) 1 [°° W+(t) 



W)= Y V +~ / d K 2 » ( 3 ) 

q 2 - m B* * J(m B +m*y t-q 2 -ie 
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where tub* = 5.325 GeV is the mass of the B*(l~) meson which induces a pole below the 
lowest multiparticle threshold at q 2 = (tub + m 7r ) 2 . Vector meson dominance assumes that 
the form factor is dominated, for all physical q 2 , by the first term in (jSJ): 

/ + (0) is the only free parameter of the VMD parametrisation. Becirevic and Kaidalov 
suggested as an alternative parametrisation the replacement of the second term in © by an 
effective pole at higher mass, 

/ ( g 2)| = ri | r2 = /+(°) | rq 2 /m 2 B , 

+ Ibz 1 — g 2 /m|» 1 — aq 2 /m 2 B * l — q 2 /m 2 B * (1 — q 2 /m 2 B *)(l — a q 2 /m 2 B ,) 

with the three parameters (ri,r 2 ,a) or (/ + (0), r, a), which are related by 

/ + (0)=n+r 2 , r = r 2 {a - 1), (6) 

where < a < 1 parameterises the position of the effective pole. This parametrisation was 
used by BZ to describe the results from LCSRs [2], whereas Becirevic and Kaidalov, faced 
with the challenge to fit three independent parameters to lattice data with limited accuracy, 
implemented the additional constraint r = af + (0) motivated from heavy quark expansion, 
and obtained the following expression in terms of two parameters, (c#, a) or (/ + (0), a): 

f ( 2xi = c B (l - a) = /+(0) (7] 

J + iy;,BK (l-q 2 /m 2 Bt )(l-aq 2 /m 2 B *) ~ (I - q 2 /m 2 B *)(l - aq 2 /m 2 B ,) ' K ) 

The BaBar collaboration has measured the integrated g 2 -spectrum of B — > niu in five 
bins in q 2 jTj, which allows one to confront the above parametrisations with experiment. 
Fitting the data by the VMD formula (J1J), we find x 2 — 9.2/4 d.o.f. The BK parametrisation 
(J7| fits the data with \ 2 = 3.5/3 d.o.f. and a = 0.61 ± 0.09. Fitting the BZ parametrisation 
(0) is slightly more subtle, as the data prefer a — » 1 or even larger, which is outside the 
allowed parameter-space. 1 Bounding a < 1, we find a minimum y 2 = 1.65/2 d.o.f. and 
r//+(0) = 0.24 ± 0.08, a = 1.00±g.is- While this implies that the VMD parametrisation is 
disfavoured, 2 both BK and BZ are viable parametrisations. Motivated by these results, we 
formulate the following strategy for extracting \V u b\ from the data: we 

• calculate /+ from LCSRs for values of q 2 where the method is applicable; 

• extrapolate the results to all q 2 using the experimentally favoured BK and BZ parametri- 
sations; 

• use the experimental information on the spectrum to constrain the input parameters 
of the LCSRs, in particular the leading-twist 7r distribution amplitude; 



1 For a = 1 the rescaling ijBJ is no longer valid. 

2 Actually VMD is disfavoured not only from the experimental point of view, but also from the theoretical 
one, as the values of /+(0) do not agree with indepedent determinations of the residue of the B* pole, see 
the discussion in Sec. 4 of Ref. 
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Figure 1: Respective size of different contributions in to f+(q 2 ) as a function of q 2 for 
central values of the LCSR input parameters. Solid line: = 0. Remaining lines from top 
to bottom (at q 2 = 0): contributions proportional to a%, dg, a%, a\. 

• determine \V u b\ from the total branching ratio using the BK and BZ parametrisations 
of/+- 

Let us start with the calculation of /+. In Ref. j2], we have presented a comprehensive 
analysis of B — > (jr, K, rj) decay form factors calculated from QCD sum rules on the light- 
cone, to 0(a s ) accuracy for twist- 2 and the dominant twist-3 contributions; earlier analyses 
can be found in Refs. jH E] . We refer to these papers for an explanation of the method. 
The main theoretical uncertainty of these analyses comes from the pion's leading-twist light- 
cone distribution amplitude (DA) <p n ; other sources of uncertainty include the b quark mass, 
the quark condensate and sum rule specific parameters (Borel parameter and continuum 
threshold). The resulting total uncertainty of /+ is between 10% and 13% j2]. Whereas the 
other parameters mainly determine the normalisation of the form factor, (p n affects also and 
in particular the g 2 -dependence and hence can be constrained from the measured spectrum. 
The DA is usually expressed in terms of its conformal expansion, 

MU, V) = 6W(1 " U) + J2 «2n(^)C2n 2 (2w - 1) ^ , (8) 

where u is the momentum fraction of the quark in the ir and runs from to 1. The C^j 2 are 
Gegenbauer polynomials and a n , the so-called Gegenbauer moments, are hadronic parameters 
which depend on the factorisation scale \i. The respective contributions of a^ <8 to /+ are 
shown in Fig. d for a typical choice of input parameters. The plot reveals that the q 2 - 
dependence of the form factor is mostly sensititive to and only to a lesser extent to higher 
Gegenbauer-moments, which agrees with the findings of Ref. ^U] . We hence decide against 
using the models for <j) n proposed in Ref. ^U], but stick with the expansion (jHJ), which we 
truncate after the contribution in a\. The values of the lowest lying Gegenbauer moments 
have been constrained from various sources, cf. Ref. [TTJ H21 El but still come with 
rather large uncertainties. A very conservative range of allowed values consistent with all 
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known constraints is 

< aj(l GeV) < 0.3, -0.15 < oJ(lGeV) < 0.15. (9) 

In this letter, we aim to constrain 0^4 from experimental data within the above range. 

We obtain values for f+(q 2 ) in dependence on a£, a\ and m& using the following criteria 
for the evaluation of the LCSRs: 

• we calculate f+{q 2 ) as a function of the Borel parameter M 2 and the continuum thresh- 
old s for two different values of m b , m b e {4.7, 4.8} GeV, five different values of 
q 2 , q 2 E {0,2.5,5,7.5, 10} GeV 2 , and 16 different values of (a£(l GeV), aj(l GeV)): 
a\ E {0,0.1,0.2,0.3}, a\ G {-0.15, -0.05, 0.05, 0.15}; 3 we interpolate the results in a\ 
and a\ in order to obtain / + as a smooth function of these parameters; 

• for each value of the input parameters, we determine /+ at the minimum in the Borel 
parameter M 2 , which implies that M 2 becomes mildly dependent on q 2 . As discussed 
in Ref. [2], this procedure ensures that a LCSR for m B , obtained from the derivative 
of the LCSR for / + in M 2 , yields the physical value m B = 5.28 GeV; 

• we choose the continuum threshold so in such a way that the continuum contribution 
to the LCSR is constant for all q 2 ; this implies that also so becomes (mildly) dependent 
on q 2 . For each value of the input parameters, we calculate /+ for three different values 
of the continuum contribution, 15%, 20% and 25%; 

• the LCSR actually yields /b/+, Jb being the leptonic decay constant of the B meson; 
in order to extract /+, we divide the LCSR by fs as calculated from a QCD sum rule 
to the same accuracy in a s . 

Some of these criteria differ from those applied in Ref. j2]. We chose to modify the criteria 
used in our previous work since we focus, in this letter, on the dependence of the form factor 
on 4, which can be meaningfully determined only if / + is calculated using exactly the same 
criteria for all values of q 2 , nib, 0-2 an d a>1- It is for this reason that we require the continuum 
contribution to be the same for all input parameters. The drawback of this procedure is that 
the sum rule specific parameters M 2 and s both become dependent on q 2 , so that, in order 
to keep the calculational effort at a manageable level, we have to restrict ourselves to a few 
points in q 2 . For each value of /+ we calculate the theoretical uncertainty by varying 

• M 2 by 40% around the central value; 

• s by ±1 GeV 2 ; 

• the central value 20% of the continuum contribution between 15% and 25%; 

• the central value of the quark condensate, — (55) (1 GeV) = (0.24 GeV) 3 , between 
((0.24 ±0.01) GeV) 3 . 

3 In the actual calculation, aj 4 are scaled up to the factorisation scale [i — \J nig — ml ~ 2.2 GeV using 
NLO evaluation. 
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Figure 2: /+ ee (0)//+ ull (0) as a function of the Borel parameter M 2 for a?, = 0.1, a\ = 0, 
m b = 4.8 GeV, s* rco = 35.25 GeV 2 , 4 ul1 = 34.42 GeV 2 , which corresponds to 20% continuum 
contribution. The ratio of the decay constants is /Ij 111 //^ 66 = 1-26. 



m b [GeV] 


4.7 4.8 


f B [MeV] 


191 ±6 164 ±4 



Table 1: fs from a QCD sum rule to 0(a s ) accuracy; the error is obtained by varying 
(M 2 )f B by 40% and (s )/ B by ±1 GeV 2 . The current world average from lattice calculations 
is f B = (189 ±27) MeV [Hj. 

The above ranges of sum rule parameters are rather conservative and account for the "sys- 
tematic" uncertainty of QCD sum rule calculations. All errors are added linearly, which 
yields a typical theoretical uncertainty of 8%. At this point we would like to comment on 
the treatment of The motivation for calculating from a QCD sum rule instead of 
using, for instance, the current world average from lattice calculations [T4] . is that (a) f B 
receives large 0(a s ) and 0{a 2 s ) corrections from gluon-exchange diagrams that also enter the 
LCSR for fsf+ and (b) fs is very sensitive to m&. By dividing the LCSR for fsf+ by fs 
obtained from a QCD sum rule to the same accuracy in a s , and using the same value of m&, 
one expects those large contributions to cancel and to reduce the sensitivity to the value of 
m&. One can check the extent to which the cancellation takes place by comparing the results 
of our 0(a s ) calculation with that at tree-level. We calculate /| ree (0) = {f B f+{0)) trce /{fB) tree 
for rrib = 4.8 GeV, a\ = 0.1, a\ = using the same criteria as for the full form factor /^ u (0) 
including 0(a s ) corrections. In Fig. |2Jwe plot the ratio /+ ce (0)//^l lll (0) for the respective 
optimum values function of M 2 . The minima in M 2 are around 7 GeV 2 . The ra- 

tio is nearly constant ~ 1.05, whereas the ratio of the decay constants, /b 111 //^ 60 is 1-26. 
This means that there is indeed a strong cancellation between the radiative corrections to the 
LCSR for /#/+(0) and the QCD sum rule for fs- We have checked that similar cancellations 
also occur for other values of mj and nonzero q 2 . Nonetheless there is a residual uncertainty 
due to the treatment of fs which we estimate to be about half the difference between the 
value of J'b calculated from a QCD sum rule to 0(a s ) accuracy and the central lattice value, 




Figure 3: Left: f+{q 2 ) as a function of q 2 for m& = 4.8 GeV, a\ = 0.115 and a\ = —0.015 
fitted to the BK parametrisation (JZj) (solid line) and to the BZ parametrisation (jHJ) (dot- 
dashed line). Right: corresponding integrated spectra in the five g 2 -bins used by BaBar. 



i.e. about 5%. The QCD sum rule results for fs are given, with errors, in Tab. ^ Adding 
these two errors linealy, we obtain an uncertainty of Jb of 8% which is independent of q 2 
and translates into a 8% uncertainty of the normalisation of / + , which we treat separately 
from the error of /#/+ calculated as described above. 4 

The next task is to compare the form factor predictions to data and to determine best-fit 
values of a£ 4 , using the experimentally favoured BK and BZ parametrisations in order to 
extrapolate the LCSR results to all physical q 2 . It turns out that the LCSR results are 
actually described extremely well by the BK and BZ parametrisations, to within better 
than 0.5%, as already noted in Ref. [2]. In Fig. |3] we show the difference between the 
BK and the BZ fit, both for the form factor and the integrated g 2 -spectra in the ith bin, 
5(i — l)GeV 2 < q 2 < 5iGeV 2 , 1 < i < 5. The form factors start to noticeably deviate only 
for very large q 2 . For the fit of the experimental spectrum, we treat the experimental errors 
as uncorrelated, but allow for a correlation of theory errors and perform the least-x 2 fits 
using the following \ 2 function and error matrix E: 

X 2 = J2(Bt h ~ BDiE-'UB? - £f p ), 

En = « XP )% + (ASf )% + C 2 {ABf){ABf )(1 - <%) , 

where Bi is the partial branching fractions ith g 2 -bin, ABf 1 the corresponding theory error 
(without the error of the overall normalisation of /+) and C is the correlation of theory errors 
which we vary within < C < 1. 

We first study the BK parametrisation which features one parameter that can be de- 
termined from the experimental spectrum: a = 0.61 ± 0.09. As / + depends on actually 
three parameters, a^, a\ and rrib, only one of them can be constrained. In Tab. |2]we give 
the best-fit values of a\ for a\ and rrib fixed. The table reveals that it is indeed possible to 
reproduce the experimental central value of a for any given a\ and rrib- It is also obvious 

4 As the quark condensate is a common input parameter for both fsf+ and Jb, the effect of its variation 
is included only once, in the uncertainty of 



6 



1. 2. 3. 4. 5. 



1. 2. 3. 4. 5. 



Figure 4: (Colour online) Best fit results for the partially integrated differential decay rates 
(solid lines) and experimental data (partially integrated spectra in g 2 -bins 1 to 5). Left: 
BK parametrisation of the form factor, right: BZ parametrisation. Solid lines: set of input 
parameters as in Tabs. |U and El Dotted line: VMD parametrisation, Eq. (0J). Dash-dotted 
line: best BZ fit to experimental data. 

that the impact of the precise value of a\ is much smaller than that of m^. Averaging over 
a\ within [—0.15,0.15] and over m?,, we find 

4(lGeV)| BK = 0.14^. (10) 

We would like to stress that this is a completely new determination of a\ and agrees very well 
with other determinations of this parameter [TTJ^JEI|. What is truly remarkable, however, 
is that, despite different best-fit values of aj, the resulting values of /+(0), and hence \V U \\^ 
agree within 3%. That is: the theory error due to and aj 4 gets largely diminished by 
the constraints on the spectrum. Using the parameter sets in Tab. El we plot, in Fig. 01 
the partially integrated spectra, normalised to the full branching ratio, together with the 
experimental data. All six parameter sets produce nearly the same curve which coincides 
with the best experimental fit using the BK parametrisation. Using the average value of the 
branching ratio as given by HFAG [IH], B(B° — > tt~£ + U£) = (1.36 ±0.11) x 10~ 4 , we obtain 

|K&|bk = (3.2 ± 0.1 ± 0.1 ± 0.3) x 10" 3 , (11) 

where the first error is experimental, the second comes from the uncertainty in the shape, 
due to the spread of values of 4 , and the third is from the normalisation of the form factor. 

Let us now turn to the BZ parametrisation, which featurs two parameters that can be 
determined from the shape of the spectrum, r// + (0) = 0.24 ±0.08 and a = 1.00±q 15 , which 
allows one to constrain for instance a\ and a\ in dependence on nib- The resulting constraints 
are shown in Fig.|51 The minimum \ 2 f° r ( a 2> a X) within the range specified in © is reached 
for = 0.23, a\ = 0.15, i.e. at the border of the parameter space; fixing a\ = 0.15, the 
best-fit value of a\ is 

4(1 GeV)| BZ = 0.23 ±0.15 for af(l GeV) = 0.15. (12) 

The contours shown in the figure include all (aj , aj) for which the fit of the corresponding 
form factor to the data yields x 2 ^ Xmm + 1 = 4.06. We immediately read off the following 
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nib 


al 


al 

A 


/+(0) 


a 


\v ub \ 




4.8 GeV 


-0.15 


0.19 ±0.13 


0.277 


0.61 


3.17 x 10" 


-3 







0.18 ±0.14 


0.272 


0.61 


3.23 x 10' 


-3 




0.15 


0.16 ±0.15 


0.268 


0.61 


3.29 x 10" 


-3 


4.7 GeV 


-0.15 


0.08 ±0.13 


0.274 


0.61 


3.21 x 10- 


-3 







0.11 ±0.14 


0.272 


0.61 


3.23 x 10- 


-3 




0.15 


0.14 ±0.16 


0.271 


0.61 


3.25 x 10- 


-3 



Table 2: 0^(1 GeV) from the fit to the shape of the experimental g 2 -spectrum, for fixed 
rrib and aj using the BK parametrisation (JZj). The table also gives the corresponding BK 
parameters /+(0) and a and the resulting central value of \V u b\. Xmm = 3.5/3 d.o.f. 





a\ 


al 


/+(0) 


r 


a 


x 2 


\v ub \ 




4.8 GeV 





0.15 


0.268 


0.18 


0.54 


3.67 


3.26 x 10- 


-3 




0.15 


0.23 


0.278 


0.14 


0.73 


3.06 


3.22 x 10- 


-3 


4.7 GeV 





0.10 


0.270 


0.18 


0.52 


3.71 


3.26 x 10" 


-3 




0.15 


0.24 


0.282 


0.14 


0.70 


3.11 


3.23 x 10- 


-3 



Table 3: Ditto for the BZ parametrisation (jSJ. In addition we give the minimum x 2 , for 
3 d.o.f., for the fit of the form factor to the data. 

constraints: 

aj(l GeV) > 0, <(1 GeV) > -0.07. (13) 

In Tab. |SJ we give the fit results for = and 0.15. The best-fit parameters do not 
agree, for within the range specified in Q, with those favoured by experiment. We have 
already pointed out earlier that the experimentally favoured value a = 1 is actually outside 
the theoretically allowed parameter space and is not supported by the results of the LCSR 
calculation. As before, we find that the different values for 034 result in nearly the same 
values of \V u b\. In Fig. 0] we plot the best-fit values for the partially integrated branching 
fractions, obtained from the parameter sets in Tab. El in comparison with the experimental 
results. For \V u b\, we find 

|K&|bz = (3-2 ± 0.1 ± 0.1 ± 0.3) x 10~ 3 (14) 

which agrees with (JTTJ). 

Combining the results from the BK and the BZ analyses, we get the following final result 
for \V ub \: 

\V ub \ = (3.2 ± 0.1 ± 0.1 ± 0.3) x 10~ 3 ; (15) 

the first error in \V u b\ is experimental, the second comes from the shape of the form factor 
and the third from the overall normalisation of / + . As for the constraints on a\±, our final 
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Figure 5: (Colour online) Allowed values of 0^4 for m& = 4.7 GeV (left) and m& = 4.8 GeV 
(right). The minimum x 2 is 3.06 and the contours include (af , aj) with x 2 < 4.06. 

results are 

4(1 GeV) = 0.19 ±0.19, aj(l GeV) > -0.07. (16) 

The value of (f 2 is the weighted average of (jlOj) and (|12|) . We would like to stress again that 
these values result from a new and independent determination of these parameters, which is 
both consistent with and complementary to the results found in Refs. [TTl 1121 ITS] . 

To summarize, we have discussed the constraints posed by recent experimental data 
on the shape of the B — > 71 decay form factor / + . We have found that both the BK 
parametrisation Eq. (JJJ) and the BZ parametrisation Eq. (J^j) can describe the data, whereas 
vector meson dominance is disfavoured. We have calculated /+ from QCD sum rules on 
the light-cone for small to moderate q 2 and extrapolated the form factor to all physical q 2 
using the BK and BZ parametrisations, respectively. We have then used the experimental 
data to constrain the pion distribution amplitude 4> n that enters the calculation of /+. We 
found that, although the best-fit values of the Gegenbauer moments 0^4 depend on mb, 
the resulting predictions for the form factors and the partial branching fractions are largely 
independent of the input parameters. This is one of the main results of this letter: the 
experimental information on the shape of the spectrum reduces the theoretical uncertainty 
of the prediction for the total branching ratio and hence the extracted value of \Vub\. In order 
to constrain <p n even further, it will be necessary to perform a combined analysis including 
also other experimental constraints from e.g. the 7r electromagnetic form factor and the 
7T-7 transition form factor [13J. 

The determination of \V u b\ presented in this letter can be improved in the future from 
both the experimental and the theoretical side. As for the latter, we would like to stress 
that the normalisation of /+ depends on the treatment of Jb, the decay constant of the B 
meson. We have determined fs from a QCD sum rule to 0(a s ) accuracy, using the same 
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values for quark mass and quark condensate as in the LCSR for /#/+, and we have shown, 
by comparison with the corresponding tree- level sum rules, that the individually large 0(a s ) 
corrections to fsf+ and js cancel in the ratio to a large extent. We have estimated the 
uncertainty of this procedure to be ~ 8%, which enters the normalisation of the form factor. 
This overall uncertainty can be reduced by calculating e.g. 0(Poa 2 ) corrections to the LCSRs, 
which is a formidable, but not impossible task. The full 0(a 2 ) corrections to fg are known 
and actually also reduce the uncertainty coming from m& [TH]- As for the experimental input, 
smaller bins in q 2 would help to further constrain the shape of the form factor and ultimately 
avoid the necessity for extrapolation in q 2 . 
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